Extended Cahill-Glauber formalism for finite-dimensional spaces: I. Fundamentals 
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The Cahill-Glauber approach for quantum mechanics on phase-space is extended to the finite dimen- 
sional case through the use of discrete coherent states. All properties and features of the continuous 
formalism are appropriately generalized. The continuum results are promptly recovered as a limiting 
case. The Jacobi Theta functions are shown to have a prominent role in the context. 
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The search for discrete quantum phase-space quasiprobability distribution functions is a subject of continuous 
and growing interest in the literature [IIIIIIIIIIESISIIIÍSIÍIIÍIEI. The possibility of representing 
quantum systems characterized by a finite-dimensional state space by such discrete quasidistributions lays the ground 
, for interesting developments and fruitful applications on quantum computation and quantum information theory 
[Ï2l IT^l IbH ITfil IÏ71 ITsl lïflj . It is well known that, as a well established counterpart to the discrete case, a huge 
variety of quasiprobabilty distribution functions can be defined upon continuous phase-space |2Cj . In this sense, the 
Cahill-Glauber (CG) approach pïj to the subject has proved to be a powerful mapping technique that provides a 
1 \ general class of quasiprobability distribution functions, where the Wigner, Glauber-Sudarshan and Husimi functions 
appear as particular cases. Therefore, it might be considered as a wide-range phase-space approach to quantum 
' mechanics regarding degrees of frcedom with classical counterparts. 

The aim of this paper is to present a discrete extension of the CG approach. Such extension is not obtaincd from 
H that approach but, instead, properly constructed out of the finite dimensional context. Furthermore, this ah initio 
^ , construction inherently embodies the discrete anàlogues of the desired properties of the CG formalism. In particular, 
discrete Wigner, Husimi and Glauber-Sudarshan quasiprobability distribution functions are obtained. Thus, besides 
the theoretical interest of its own, such extension has direct applications in quantum information processing, quantum 
tomography and quantum teleportation, which are explored in a following work . 

This work is organized as follows: In the next section we briefiy outline the CG approach, setting the stage for 
section III, where our proposal for a discrete extension of the CG mapping kernel is presented. In section IV bàsic 
properties of the mapping technique are discussed, and the continuum limit is carried out on section V. Finally, section 
VI contains our summary and conclusions. Also, important calculations are detailed in the Appendix. 



II. THE CAHILL-GLAUBER MAPPING KERNEL 



For the sake of clarity, in what follows we will briefiy review the central ideas which constitute the core of the CG 
approach, and that will be properly generalized in the following sections. Basically, the cornerstone of the formalism 
is the mapping kernel (hereafter h = 1) 



T W (í,p) 



dp'dq' 
2tt 



exp [— ip'{q - Q)] exp [iq'(p — P)] exp ( —-p'q' ) exp -(ç' 2 + p' 2 ) 



(D 



which is responsible for the mapping of bounded operators on the continuous phase-space, being s a complex variable 
satisfying the condition \s\ < 1. Here, the momentum and coordinate operators obey the Weyl-Heisenberg commuta- 
tion relation [Q, P] = il. Since the above expression explicitly depends on s, this parameter labels an infmite family 
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of mapping kernels. Each mapping kernel can be seen as the double Fourier transform of the displacement generators 
multiplied by a phase factor exp [(i/2)p'q'} and by the folding function exp [(s/4)(q' 2 +p' 2 )] . For purposes which will 
become evident later, we write the mapping kernel as 

T W (ï,p) = J°° ^^exp[-V(?-Q)]exp[V(p-P)]exp(-Í p y) ((0\q',p'))- s » (2) 

where \q',p'} is a coherent state. 

The mapping of a given operator is achieved by the trace operation 0^(q,p) = Tr[T^ s \q,p)0], being <D^(q,p) 
the function which represents O in the associated usual phase-space. The mapping is one-to-one, and the operator is 
reobtained from its associated function by 

O = /"^0W( ít p)T<-)( g> p). 

It is clear that, for each operator, there is an infinite family of associated functions labeled by s. In particular, the 
phase-space representatives of the density operator are referred to as quasiprobability distributions functions and 
have, obviously, distinguishable importance [2(J. One of the great virtues of the Cahill-Glauber approach is that three 
special and important types of quasidistributions, namely the Glauber-Sudarshan (s = 1), Wigner (s = 0) and Husimi 
functions (s = —1), are particular cases. Each of these functions have been extensively explored and reviewed in the 
literature |23 |. 

A particular mapping kernel, characterized by a given parameter s, can be expressed in terms of another mapping 
kernel with a different parameter value. The same holds true to the functions associated with a given operator. In 
fact, the procedure in the latter case can be easily shown to be the same as in the former. That is, we may discuss only 
the relation between the mapping kernels, knowing that equivalent relations are observed by the associated functions. 
In this way, the connection between the two mapping kernels is seen to be given by the trace of the product 



Tr[T^)( gi , Pl )T^)(g 2 ,p 2 )] = ^ ^exp{i[q(p 1 -p 2 )-p(g 1 -g 2 )]}exp 



«o|?,p»- (si+S2) 

(3) 



— : — (<? +P ) 



Sl + S-2 



exp(— [(p! - P2 ) 2 + (qi -ç 2 ) 2 ] j Re( Sl + s 2 ) < . (4) 
l si + s 2 L J 



We immediately recognize the important role played by the last exponential function in J3J), since, if the condition 
Re(si + s 2 ) < is not observed, the trace gives a divergent result. Thus, that condition imposes a constraint that 
defines a hierarchy. That is, on continuous phase space there is a hierarchical structure of mapping kernels allowing 
one to express a given phase-space function in terms of a Gaussian smoothing of another, and, as such, inverse relations 
do not exist. In other words, the Gaussian folding hierarchical structure observed by the quasidistributions has its 
roots in the functional form of (0\q,p). We stress this particular point as the discrete equivalent to equation @ does 
not imply a hierarchical relation. 



III. THE DISCRETE MAPPING KERNEL 



A. Preliminaries 

1. Operator bases 

Long ago Schwinger proposed the following set of operators to act as a basis on an operator space 

s(^) = ^un^cx P (^) , 

where the U's and V's are the so-called Schwinger unitary operators [24[, N is the dimension of the associated state 
space and the indices {??,£} run on any complete set of residues mod(iV); in particular we choose the closed interval 
[—£,£], with i = (N — l)/2. For simplicity, we shall restrict ourselves to the odd N case. Even dimensionalities, for 
the purposes of this paper, can also be dealt with simply by working on non-symmetrized intervals. 
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The set {S(rj, Ç)} v spans a complete and orthonormal basis on the N 2 space of linear operators acting 
on finitc complex vectorial spaces, in the sense that, as the trace operation stands as the inner product on operator 
spaces, any linear operator can be written as 

t 

o= r &[stfa,Oo]sfa > 0· (5) 

The fundamental result 

ensures that this decomposition is unique. The superscript [N] on the Kroenecker deltas denotes that they are different 
from zero whenever their indices are mod(-ZV) congruent. The Schwinger basis elements also obey the property @ 

S+(7 ? ,£) = S(-7 7 ,-0· (6) 
2. Discrete Coherent States 

The Schwinger operator bases elements also act as displacement operators on a particular reference state to form 
discrete coherent states as 



|r ? ,0 = ViVS( ? 7,-0|0,0), (7) 

where the reference state is written by means of the Jacobi í? 3 -function (whose explicit form is shown in Appendix 
A) as 

1 1 

l°'°) = X H #3 (2a 7 |2ia) K) , (8) 

7=-* 

where {|u 7 )} 7= _^.. are the eigenstates of the unitary operator U, 

= -1= [i? 3 (0|ta)i?3(0|4io)+i?4(0|ta)i?2(0|4io)] 

is the normalization constant, and a = (2N)^ 1 . Due to the properties of the i?3-function, the reference state above is 
preserved under the action of the Fourier operator 0, [2í| 



310,0) = |o,o) , 

where 

t 

l U 7>(%l » 

and {|w 7 )} 7 =-£,...,f are the eigenstates of V, with (u^v^) — exp[(2Tri/N)(jrf\. Parity of the í?3-function also ensures 
that (u K \Q, 0) = (u- K \0, 0), from which it follows 

(0,0|/^> = <0,0|-M,-^>. (9) 

There are, of course, a number different recipes of discrete coherent states, some of them also in connection with 
i9-functions, for instance [26ll27|| . 

B. The extended mapping kernel 

Now let us define the extended mapping kernel as 
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where s is a complex number satisfying \s\ < 1, and %(rj, £) = (0, 0\r], £) denotes the overlap of coherent states explicitly 
calculated in Appendix A. The set {S(?7, C)}n,í=-i ... £ itself spans a complete and orthogonal basis on operator space. 
Nevertheless, we can go back to decomposition {JJ, use equation ©, and introduce convenient factors to get 

l 

0= Yl Tr [S(-7?, -flO] Sfa e) [JC(-7?, -Q] s [Xfa 0]' s , 

rç,£=-« 4 : ' 



where equation JÏÏJ has been used. Conveniently grouping the terms the new decomposition reads 

t 

°= E ^ [ s(_s) Hj.-Oo] s« ■ 



(10) 



Now, introducing the double Fourier transform of (77, £), i.e 

e 

i 



tW(t 7> = -^ E S(s) (^6exp 

fJ,,l/=—& 

and its Fourier inverse, we can, after a few steps, write equation I|1U|) as 



(11) 



with 0^ 3 >(ii,v) — Tr [T^ s )(/í, ;/)0], defining a one-to-one mapping between operators and functions defined on a 
discrete phase-space {/i, v}, where explicitly 



2-kí 



N 



exp(-^) [X(V,0]~ S , 



tU(jí,v) = - J2 U"V«exp 

n>í=-£ 

and 3C(r), £) can be shown to be a sum of products of Jacobi t9-functions (as seen in Appendix A), 

1 



(12) 



{ r &3(ar]\ia)i!)3(a£ l \ia) + ·d3(ar)\ia)'&4,(aÇ\ia) exp(z7ríy) 



+ í?4(a?7|m)í?3(aÇ|ia) exp(Í7r£) + i?4(aí7|ia)i?4(fl£|ifl) exp [ÍTr(rj + £ + N)]} 



(13) 



The new kernel, written as in equation l|12(l . allows us to conclude that the above sum of products of ^-functions 
plays, in the discrete phase-space, the role reserved to the Gaussians in the continuous case. 



IV. PROPERTIES 



A. Basic general properties 

From the properties of the mapping kernel it is straightforward to obtain general properties of the associated 
functions in phase-space. We observe that all following properties correctly generalize the continuous CG ones. First 
we note that 



(i) T«( M ,z/) =T( S *)( M!Í ,) 



(14) 



implying that the mapping kernel is Hcrmitian for real vàlues of the parameter s. As a direct consequence, the 
phase-space representatives of Hermitian operators are real. 
Direct calculations also show that 



(ü) l Yl t(í) (^) = i. 



= i 



(iii) Tr[T^ )V ) 

(iv) Tr TUfoi/jTt-'Vy)! =Nól N l4 N l j, 



(15) 

(16) 
(17) 
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The property (iv) is a crucial one from which expression Qllj l could be immediately obtained. From tiris property also 
follows the general result 



Tr(AB) = Í fi ^"W)^ -0 (*»>")■ 

fi,v=-í 

In fact, property (iv) is a particular case of the general expression 

1 



Tr 



tM( M!Í /)TM( M >') 



N 



E 



»7,í=- 



(18) 



which is the counterpart of equation |J3J. It must be stressed that this expression is always well defined, even for 
Re(í + s) < 0, as X(rj, £) ^ 0. 



B. Particular cases 



There are three important particular cases to be discussed: 
i) (s=0) In such a case it is easy to see that 



(19) 



where G(/i, v) is the mapping kernel introduced by Galetti and Piza, which is a discrete gencralization of the 
Weyl-Wigner mapping kernel 0,0, EEH^!- in that case, being p the density operator, equation (fTT|> would 
read, for O = p, 



H,v=—t 

with W(/z, ^) = Tr [G(/i, v)p\ a discrete Wigner function. 
ii) (s=-l) A fundamental property of our mapping kernel is 



(20) 



which can be proved decomposing the coherent state projector in the Schwinger operator basis as (using equations 
© and ©) 



M ,z/)( / i,HV- í U-"exp(--ry£ 



\V,u)(iM,v\ = jï E U"V«expfeATr 

which, by applying the defmition of the coherent states, equation {JJ, reads 

1 £ 

\w)^v\ = ü E u"v«(o,o|v-u^v-íu-"U"v-ho,o), 



»7,í=- 



and using the Weyl commutation relation, U a V^ = exp[— (27Tï/A f )a/?]V' 3 U c 



.">M = jy E U"V«exp 



i7,e=-« 



2ttï , 



exp( -r?£) (0,01*7, f) 



where in the last step parity of (0, 0|?7, £) with respect to 7/ was used. This proves our assertion. As a consequence, 
the phase-space decomposition of the density operator, associated with this particular value of the parameter s, 
reads 



1 

H,u=—t 

allowing us to identify P(/i, v) as a discrete Glauber-Sudarshan distribution. 



iii) (s=l) In this case we may write 



P=N E ^ { - l) Mp\^ {1) {^) 



n,v=-e. 



which is simply 



H,u=—t 

By definition JC(/i, v) — (fi, v\p\[í, v) is positive definite, and it can be idcntificd as a discrete Husimi function. 
As any operator can be decomposed by the use of expression , it follows that we are allowed to write 

(7,A= — & 

where the minus signal was kept only for clarity. We then use equation l|18|) to write explicitly 



Tr 



»?,í=- 



(21) 



that is, the discrete Fourier transform of the %(r],£) is the folding function. The above result can also be written in 
the compact form 



Tr 



which is precisely the Wigner function associated with a coherent state v). We therefore have 

1 1 

T ( - 1) (M,^) = T7 £ ( m ,Hg(^a)I/^}t(°V,a). 



(22) 



er,A=-£ 



In the same form we now decompose T(°)(^t, v) as 



T(°)( M ^) = Í £ Tr TCV^T^W) T«(<t,A), 



which allows us to use once again the above result for the trace and write 

l 



T<°>(a*,i/) = Í fi (M,HG(tr,A)|M,i/)T( 1 )( £ r,A) 



(23) 



cr,A=-£ 



Multiplying both equations (|2*2*|l and (|23|l by the density operator p and taking the trace, we are led to the suggestive 
results 



1 - 

= ^ E ^,f|G(<7,A)|//, I /)W( < 7,A) , 

cr,A=-^ 
1 ' 

w (^) = n E (m^|g(^a)|m,^a), 



,A=-Í> 



which are the discrete counterparts of the well known Gaussian smoothing that occurs in the continuous case, in 
agreement with the hierarchy present in that context. 
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It must be stressed, however, that, opposed to the continuous case, it is now possible to write 

1 1 



i 



T ( %,^) = \ È A( At - < 7,i/-A)T< >(<7,A), 



-,A=-t? 



where 



»?>«=--< 



2ttí 



A(/x-<r^-A) = - ]T exp^— [r,( M -a)+^-A)]^[aC(r ? ,0]" 1 , 



TV 



(24) 



which, at least in principle, can always be calculated (we remind again that 9C(r),t;) is frnite and diffcrcnt from zero). 
Also very illustrative is the result that follows from the decomposition 



(25) 



With 



Tr 



N 



2-kí 



J2 ex P \—l V (n-<j)+av-V}\inv,0} 



T(- 1 '( ( r,A)T(- 1 )( M , í ,) 
which can be shown to be v\<j, A)| 2 , we have 

i f 

T(- 1 )( M)I /) = - £ K^M^T^A) 



(26) 



Thus t/]cr, A)| 2 itself, which is the Husimi function associated with the discrete coherent-state \fx, is), acts here as 
the smoothing function. 



V. CONTINUUM LÍMIT 

Following the procedure detailed in both [2^, [2íj , the continuum limit of the mapping kernel I|i2(l is reached as 
follows: we introduce the scaling parameter e = (2tt /TV) 1 ' 2 , which will become infinitesimal as N — > oo, and the two 
Hcrmitian operators 

i £ 
P= E / í£ Po|^)( w mI Q= E M' e ?olv)(vl > ( 27 ) 

constructed out of the projectors of the eigenstates of U and V. The parameters po and qo, with p q = H = 1, are 
chosen to be real, carrying units of momentum and position, respectively, while epo and eqo are the distance between 
successive eigenvalues of the P and Q operators. Then, rewriting the Schwinger operators as 

V = exp(^) U = e X pf^V (28) 

and performing the change of variables q' — — <Zo e £i p' — Po^li P — Po^ v and q = Ço^M? we obtain 

T W (?,P)= E E ^^expHp'(9-Q)]expKb-P)][WPoe,-979oe)]- s expí~^V 

q' = -q eep' = -p el 

As TV — > oo, it follows that Ao' — * cfg' and Ap' — ► dp'. Since the continuum limit of the discrete coherent-states has 
been already discussed in 0,E|> it is clear that the term [X(p'/poe, —q'/qoe)]~ s , which is even, will go to (((%', p'))~ s . 
Therefore we end up with 

T«(g,p) = |°° ^^exp[-ip'(q-Q)}exp[iq'(p-P)]e X p(- l -p'q^ «0|^,p'»— , 

which is exactly the mapping kernel @ of Cahill and Glauber. 
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VI. CONCLUDING REMARKS 



The results obtained here show a genuine discrete mathcmatical structure which closely parallels the one of Cahill 
and Glauber. This was achieved pursuing the lines proposed in |2j, which makes use of the discrete Fourier transform 
of the Schwinger operator basis, to deal with the discrete phase-space problem. Now, we stress that expression J5j is 
as simple as it is important, since it clarifies the role of the coherent states overlap within the CG approach. By its 
turn, the discrete coherent states proposed in [3, 1(1 l2tl |27| provide a natural path for a discrete extension of the CG 
formalism, while the properties of these states have played a crucial role as they led, for example to the bàsic equation 
(|2Up. Thus, the coherent states overlap can be seen as the link between the discrete and continuous approaches. 
Furthcrmore, the continuum limit presented in section V ensures that the CG mapping scheme is correctly recovered 
through a limiting procedure which is mathematically consistent j2^, El Ell H3| • 

It is worth mentioning that Opatrny et al Q have pursued a goal similar to ours. Although both approaches share 
virtues, our formalism presents mathematical features that allow us to achieve farther reaching results. It is precisely 
the correct choice for the reference state, and the mathematical procedure adopted here, that lead to the obtention 
of such a wide set of important results. 

The use of Schwinger operators is crucial if one is concerned with the problem of ordering. As they are unitary shift 
operators, equation (|12[1 makes it clear that the associated expansion is necessarily linked to a particular ordering of 
U and V, which can be directly connected to the Q and P ordering of the continuous case. 

Concerning the role of the Jacobi Theta functions in the discrete phase-space context - they are implicit in the 
3C(?7, £) term -, comparison with the usual CG results makes it evident that the Gaussian (or anti-Gaussian) terms, 
which are present in the continuous case, are here replaced by the sum of products of functions l|13|) . and its Fourier 
transform l|21[l ; both play here the role of the smoothing functions. 

It is always important to emphasize the discrete case's peculiar features that do not have correspondence in the 
continuum. A plain example of these is expressed by the well-behaved function given by equation 1)24(1 ■ whose 
continuum limit clear ly diverges, as the hierarchical structure presented in |0J would imply. The finite character of 
the discrete scenario prevents such a behaviour since, even if some terms in equation (|24|l might become large for large 
N, they remain always finite due to the behaviour of the z9-functions. This allows one - to give a extreme example - 
to express, in the discrete scenario, the Glauber-Sudarshan function in terms of the Husimi function. 

Finally, it is worth mentioning that the mathematical formalism developed here opens new possibilities of inves- 
tigations in quantum tomography and quantum teleportation. These considerations are under current research and 
will be published elsewhere (for instance, see reference p2^~). 
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The vacuum associated with the discrete coherent-states {7J , given by equation (jHJl , is written in terms of the Jacobi 
z?3-function, which reads 
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APPENDIX A: THE DISCRETE COHERENT-STATE OVERLAP 




with a = (2N) 



Explicitly, then, the overlap (rj, Ç\fj,,v} is: 




tf 3 (2H2m)tf 3 (2H2m)(tv|U^Wí-> K ) , 



|0,0) 



K,V = -Í 
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which, by evaluating the matrix element within the summation, gives 
1 1 

~2 exp {2iria[— \w + rjÇ + 2£(/i — rj)]} d 3 {2an\2ia) , d 3 {2a{K + n ~ rf)\2ia) exp [—irianiy — £)] 



N 2 

Now let us call 



í/)= i?3(2a/í|2ia)i?3(2a(/í + /z)|2ia) exp(— -nianv) , (Al) 



so we write 



(«7, Cl/", v) = ■L· exp {27ria[-^í/ + r)£ + 2£(/z - t))]} A(fi - rj, v - £) . 

The task now is to evaluate explicitly the term A(p, v). Prior to that, some properties of A(fi, v) might be verified 
from the start. First, we note that A{p,,v) = A(v, fi) (using the fact that the $3 is an eigenfunction of the discrete 
Fourier transform), and also A(fi + Nk, v + Nu') — A(/i, v\ for k and k' integers. 
Let us now evaluate A(fi, u). We start using the $3 deíinition to get 

í 00 

A(/i, v) — exp(— 4irianfi) exp(— 2naa 2 + Airiana) exp [— 27ra/3 2 + 4iria(K + v)0\ ■ 

K— — Í a,/5— — 00 

The sum over k can be readily carried out, which gives us NS^j l _ a . Thus, (3 will assume the vàlues /i — a + -/V7, 
where 7 is an arbitrary integer, yielding 

00 

N exp(— 2naa 2 ) exp [— 2na(pL — a + -/V7) 2 + kïïiav{p, — a + iV-y)] . 

a.7— — 00 

The last term in the second exponential is equal to one, and we shift the sum over a by /1 in order to get 



N ^2 ex P 

Q,7 — — OO 



-4ira ( a — ) 2-Kaji — A-KÍaa(v — i/i) 



Now, we split the sum over 7 in contributions coming from the even (e) and odd (o) integers as follows: A = A e + A . 
Consequently, the even term can be dealt with by shifting the sum over a by — ÍV7, 

00 

A e (/i,z/) = N exp [-Airaa 2 - 2-k^ 2 N - 2nafi 2 + íniaÇa + N^){v - 

a,7 — — cjo 

oc 00 

= N exp [— Airaa 2 + Anidaiv — i/i)] exp [— 2ir r y 2 N + &mar{N{y — ifijj 

a— — 00 7 — — oo 

= iVí?3(2a(^-^)|4m)í?3(i/i|2iiV)exp(-27ra/i 2 ) , 
and, in a similar fashion, the odd term gives 

A (p, v) = N-dz(2a[v - í(jjl + N)] \4ia) d 2 {i^\2iN) exp f _ -y- - 7r/i - 2nafi 2 

Next, using a relation (all the following relations for í?-functions come from |3^|) 



T^j = exp (-i 7r - - iirçj #2(ç|r) 



we get 

A a (fi, v) = Ntf 3 (2a[is - + N)]\4da) & 3 (i(fJ, + N)\2iN) exp [-2wa(/j, + N) 2 ] 
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Thcn 



A(p,v) = Nd 3 (2a[v - ii4\4ía) 3 (ifj,\2iN) exp (-2a7r^ 2 ) 

+ Né 3 (2a[v - i{n + N)]\Aia) i? 3 (i(/x + N)\UN) exp [-2na(fj, + N) 2 ] 

Now we recali the fundamental property of $-functions, 



^3 I - 



T \ T 



.-2 



Vrexp ( — ) i? 3 (ç| ir) , 



IT 

which applicd to the d 3 (í^\2iN) and # 3 (i(/z + iV)|2iAT) terms above gives 

j4(/x, i/) = \j^-'d 3 {2a{v — ifi)\Aia) ■d 3 (afi\ia) exp (— na^ 2 ) 

— üz{2a[v - i(fi + N)]\4ia) 3 (o(a» + W)|ia) exp [-7ra(/i + N) 2 ] . 
Still, onc has to use the relation 

*« t >=sKíID + *(IIï)] 

in the first i?-function of each of the two terms above to finally employ the relations of quasi-periodicity 

í9 3 (ç + mT|r) = exp (— nÍTm 2 + 2-ïïimç) i? 3 (ç|r) 

^4(ç + mr|r) = exp (— -Kirm 2 + 2irim<; + irim) f?4(ç|r) , 

to reach the desired result 

A(fi, v) = _ exp(27ria/w) {í? 3 (a/i|ia) , i? 3 (a^|ia) + i?3(<z//|ia)i?4(aï/|ia) exp(Í7r/i) 
4y a 

+ t?4(a/i|ia)^ 3 (o^|zo) exp(iivv) + i?4(a/i|ia)í?4(a^|ia) exp[í7r(f + fi + N)]} . 

This completes the calculation of (rj, £|/z, f). 

Finally, it is convenient to write down the result for the particular case 9C(fj,, v) = (0, 0|/U, f), i.e., 

9C(/z, i/) = ^Ij exp (—2ma/j,v) A{fi, v) . (A2) 
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